function hs = hrmsum(N)
hs = sum(1./(1:N));
end